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The scaling properties of elliptic flow of hadrons produced in ultrarelativistic heavy-ion collisions 
are investigated at low transverse momenta, pr < 2 GeV. Utilizing empirical parameterizations of a 
thermalized fireball with collective-flow fields, Resonance Recombination Model (RRM) is employed 
to describe hadronization via quark coalescence at the hadronization transition. We reconfirm that 
RRM converts equilibrium quark distribution functions into equilibrated hadron spectra including 
the effects of space-momentum correlations on elliptic flow. This provides the basis for a controlled 
extraction of quark distributions of the bulk matter at hadronization from spectra of multi-strange 
hadrons which are believed to decouple close to the critical temperature. The resulting elliptic flow 
from empirical fits at RHIC exhibits transverse kinetic-energy and valence-quark scaling. Utilizing 
the well-established concept of sequential freeze-out we find that the scaling at low momenta can be 
extended to bulk hadrons (ir, K, p) at thermal freeze-out, and thus result in an overall description 
compatible with both equilibrium hydrodynamics and quark recombination. 
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I. INTRODUCTION 

Ultrarelativistic heavy-ion collisions (URHICs) enable 
the creation and study of superdense strongly interacting 
matter, possibly associated with the formation of novel 
phases where partons deconfine and chiral symmetry is 
restored [l|. One major finding by the experimental pro- 
gram at the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory is the large azimuthal 
anisotropy of hadron transverse- momentum (pt) spec- 
tra in non-central nuclear collisions, the so-called elliptic 
flow [2j. It is quantified by the second harmonic coeffi- 
cient, V2(pt), in the Fourier expansion of the azimuthal 
angle, <f>, for the pr-spectra of hadrons. 

Experimental measurements suggest that the behav- 
ior of v-2{pt) may be classified into three regimes. At 
high transverse momentum, px > 6 GeV, the azimuthal 
anisotropy is believed to be due to the path-length de- 
pendence of high-energy partons as they traverse the hot 
medium and lose energy. We will not be concerned any 
further with this mechanism in the present work. At low 
transverse momentum, Pt ^ 2 GeV, i>2 increases with px 
essentially linearly but with a delayed onset for hadrons 
with increasing mass, giving rise to the so-called mass 
ordering of elliptic flow [3] . This regime is well described 
by hydrodynamic simulations and thus indicates a large 
degree of equilibration, encompassing more than 90% 
of the observed particles At intermediate trans- 

verse momenta, 2 GeV< pr < 6 GeV, V2(pt) saturates 
and exhibits a remarkable scaling between baryons and 
mesons j9|, 



v 2 {pt) = nqV^'ipr/riq) 



(1) 



where n q denotes the number of valence quarks contained 



in hadron h. Equation (fT]), also known as "constituent- 
quark number scaling" (CQNS), has been successfully 
described in the framework of quark coalescence mod- 
els, where quarks recombine into hadrons at the tran- 
sition from partonic to hadronic matter [Icj . In line 
with Eq. (fTJ), these models thus imply that the observed 
hadron elliptic flow can be reduced to a universal quark 
elliptic flow at hadronization. This assertion, however, 
requires several simplifying assumptions. First, inter- 
actions in the subsequent hadronic evolution (between 
hadronization at T c ~ 170 MeV and kinetic freeze-out 
at T fo ~ 110 MeV) are neglected. Furthermore, Eq. flU 
can only be dervied in instantaneous quark coalescence 
models where energy conservation is violated and narrow 
hadron wave functions in momentum space are utilized. 
In addition, the underlying quark distribution functions 
are assumed to factorize in coordinate and momentum 
space, implying that the elliptic flow of the system is 
implemented point- by-point (i.e., locally), rather than 
through a space-dependent flow field which is a conse- 
quence of a collectively expanding source. The imple- 
mentation of realistic space-momentum correlations into 
quark coalescence models has indeed been proved diffi- 
cult [ll|. Finally, quark coalescence models typically do 
not comply with the principle of detailed balance, imply- 
ing that the proper equilibrium limit cannot be estab- 
lished. This hampers attempts to extend the description 
into the 1ow-pt regime where thermalization is believed 
to prevail. This has become a more pressing issue once 
it was realized that CQNS can be generalized to encom- 
pass both low and intermediate px if «2 is plotted ver- 
sus transverse kinetic energy, KEt = ttit — m, where 
rriT = (p T + to 2 ) 1 / 2 is the hadron's transverse mass [I2T - 
IT7I ]. In this representation, elliptic- flow data for all ob- 
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served hadrons fall onto one universal curve for V2/n q 
versus KE T /n q (HHl. 

Progress on the above issues has been recently re- 
ported by reformulating the quark-coalescence process 
at hadronization in terms of resonant qq — > M scat- 
tering (M: meson) [l9j . Implementing the under- 
lying cross section into a Boltzmann transport equation 
(e.g., in Breit-Wigner approximation) not only guaran- 
tees energy-momentum conservation but also satisfies de- 
tailed balance and thus leads to the correct thermal- 
equilibrium limit. It has been verified [l8j that the me- 
son spectra resulting from the Resonance Recombination 
Model (RRM) are compatible with CQNS, Eq. (!]), in 
the factorized (local) approximation to V2(pt)- A further 
step has been taken in Ref. [2(| where the RRM has been 
applied to more realistic quark distribution functions as 
generated by relativistic Langevin simulations within an 
expanding QGP fireball [21] for semi-central Au-Au col- 
lisions at RHIC. Since the Brownian-motion approxi- 
mation underlying the Langevin process is only reliable 
for relatively massive and/or high-momentum particles, 
nix 3> T, the results have been restricted to charm 
and strange quarks. Under inclusion of the full space- 
momentum correlations imprinted on the quark distri- 
butions by the hydro-like flow fields of the expanding 
fireball, the elliptic flow of </> and J/tp mesons was found 
to exhibit CQNS in KE T from to -3 GeV, which en- 
compasses both low and intermediate pt, i.e., the quasi- 
thermal and kinetic regimes. However, several impor- 
tant questions remain, e.g., the manifestation of CQNS 
and KEt scaling for bulk (light) particles in the low-pr 
(thermal) regime, the robustness and generality of the 
flow fields utilized in the fireball simulations or the role 
of reinteractions in the hadronic phase. A thorough un- 
derstanding of these issues, in connection with a quanti- 
tative description of hadron data, is hoped to ultimately 
enable the extraction of quark distribution functions just 
prior to the conversion to hadronic degrees of freedom, 
and thus to quantitatively establish the presence of a col- 
lectively expanding partonic source in URHICs. 

In the present paper we conduct investigations of scal- 
ing properties of hadron spectra by focusing on the low- 
Pt regime coupled with the assumption of complete ki- 
netic equilibration. We will utilize the RRM to con- 
vert locally equilibrated quark distribution functions into 
hadron spectra using a blast-wave type quark source 
with realistic flow fields at the hadronization transi- 
tion (including space-momentum correlations character- 
istic of hydrodynamic simulations). Contrary to ear- 
lier work [20], we do not attempt to generate the flow 
fields from a dynamic evolution, but rather extract them 
from empirical fits to hadron p^-spectra and elliptic 
flow. In fact, after verifying that resonance recombina- 
tion converts collective, locally equilibrated quark dis- 
tribution functions into equilibrium hadron distributions 
with identical (space-momentum dependent) collective 
properties, one can adopt the source parametrization at 
the hadron level and thus readily extend the description 



to baryons. For this part of the investigation we focus on 
multi-strange hadrons (i.e., hadrons containing at least 
2 strange and/or anti-strange quarks) which are believed 
to kinetically decouple close to the hadronization tran- 
sition. This notion is supported by empirical blast-wave 
fits to RHIC and SPS data [J, as well as by 

the putative absence of resonances that these hadrons 
could form in interactions with bulk particles such as 
pions, kaons or nucleons (hadronic resonances drive the 
collective expansion of the hadronic phases in URHICs). 
We refer to hadrons decoupling close to T c as "group- 
I" particles, which are thus the prime candidates to in- 
fer properties of the quark phase just above the critical 
temperature. On the other hand, the most prominent 
bulk particles (n, K, p, etc.) are well-known to inter- 
act quasi-elastically via strong hadronic resonances, such 
as p(770), If* (892) and A(1232), which renders their in- 
teractions operative until a kinetic freeze-out tempera- 
ture, Tfo, which is significantly lower than the one of 
group-I particles, T fo ~ 110 MeV < T c ~ 170 MeV. We 
refer to hadrons decoupling at Tf Q as "group-II" parti- 
cles. Clearly, a description of KEt scaling in the \ow~pt 
regime would be incomplete (to say the least) without 
7T, K and p. In an attempt to include these particles, 
we adopt the simplest possible scenario of introducing a 
second source for group-II particles by fitting their spec- 
tra and V2 with a space-dependent elliptic flow field at 
Tf . We then evaluate in how far such a minimalistic (but 
in our opinion well-motivated) sequential freeze-out can 
be compatible with CQNS and KEt scaling, and which 
mechanisms could be responsible for them. In doing so 
we also revisit the question of resonance feed-down con- 
tributions to the pion vi for which differing results seem 
to exist in previous works. 

Our article is organized as follows. In Sec.|H]we reca- 
pitulate basic features of the RRM, specifically its equi- 
librium limit in the presence of anisotropic flow fields. In 
Scc. lIIII wc quantitatively construct realistic flow fields for 
semi-central Au-Au collisions at RHIC using group-I par- 
ticle freeze-out close to the hadronization transition. In 
Sec. IIVI we reiterate this procedure for group-II particles 
at kinetic freeze-out. In Sec.[V]we examine our results for 
group-I and -II particles with respect to CQNS and KEt 
scaling. The effects of feed-down on the inclusive pion V2 
are revisited in App. [X] We summarize and conclude in 
Sec. ED 



II. RESONANCE RECOMBINATION MODEL 
AND EQUILIBRIUM LIMIT 

Early quark coalescence models which employed in- 
stantaneous approximations in the recombination of 
quarks [ToL l24j were quite successful in providing an 
explanation for two rather unexpected phenomena in 
hadron spectra at intermediate pt in Au-Au collisions 
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at RHIC. 1 Specifically, these were the enhancement of 
baryon-to-meson ratios (B/M — p/tt, A/if) over their 
values measured in p-p collisions, as well as CQNS of 
V2(pt)- While an enhanced B/M ratio is also a nat- 
ural outcome of hydrodynamic flow, the latter is not 
expected to generate enough yield to dominate hadron 
production at pt > 2 — 3 GeV. The instantaneous pro- 
jection of parton states onto hadron states conserves 3- 
momentum by construction but energy conservation is 
violated due to the sudden approximation jTT[. This ap- 
proximation is believed to be viable at intermediate pt 
where corrections are expected to be of order 0(Q/pr) 
with Q = m — m q — niq (m: meson mass, m^n: anti- 
/quark mass), albeit with unknown coefficient [18|]. How- 
ever, this poses a significant problem at low pt and/or 
for hadrons (resonances) with large binding energy (Q- 
value) . 

In Ref. [lH quark coalescence in an interacting medium 
in the vicinity of the hadronization transition was reinter- 
preted as a process akin to the formation of resonances, 
q + q ^ M, and implemented via a Boltzmann equation 

P »dJ M (t,x,p) = -mTf M (t,x,p) +p°p(x,p), (2) 

where /m(^, x,p) denotes the phase space density of the 
meson, and the gain term is given by 



d 3 pid 3 p 2 

—fq{x,P)fq{x,p) 



(2tt) 

xa(s)v rc i(pi,p 2 )S 3 (p -pi -p 2 ) 



(3) 



The use of an explicit (resonant) cross section automat- 
ically satisfies energy-momentum conservation. For sim- 
plicity it has been modeled by a relativistic Breit-Wigner 
form, 



cr(s) = g a 



4-7T 



(Tmf 



k 2 (s — m 2 ) + (Tm) 2 



(4) 



where the same reaction rate T is used as in the loss term, 
the first term on the right-hand-side (rhs) of Eq. ©, 
and g a is the degeneracy factor. This guarantees de- 
tailed balance, which is essential for recovering thermo- 
dynamic equilibrium in the long-time limit, At 3> 1/r. 
If hadronization is rapid enough to produce hadrons in 
equilibrium, as it seems to be the case for the bulk of 
hadrons in URHICs, this limit is applicable and the re- 
sulting hadron-momentum distribution is given by 



Em(p) 
mT 



d 3 x/3(x,p) 



(5) 



For more details, we refer the reader to Ref. [18|. It 
is well known that the unique equilibrium solution of a 



1 Throughout this manuscript, we will use px (pt) to denote 
hadron (quark) transverse momenta. 
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FIG. 1: (Color online) Comparison of <^-meson pr spectra 
computed from (i) resonance recombination of strange quarks 
using blast-wave flow fields for strange quarks (rectangles), 
and (ii) a direct application of an equilibrium blast-wave dis- 
tribution (solid line) with the same flow field as for quarks 
(solid line). The hadronization temperature has been cho- 
sen at T c = 180 MeV and the radial flow at the surface as 
v = 0.55. 



Boltzmann transport equation is a Boltzmann distribu- 
tion if and only if energy conservation as well as detailed 
balance are satisfied. The RRM complies with these re- 
quirements. In the remainder of this section we will verify 
numerically that Eq. ([5]) recovers the thermal Boltzmann 
distribution for mesons in the presence of a spatially de- 
pendent anisotropic flow field at the quark level (in pre- 
vious work [3, [20| , this has only been exemplified for pr 
spectra, not for v 2 ). 

As a first test we study <f> mesons in an azimuthally 
symmetric and longitudinally boost-invariant fireball. 
The phase space-densities of the input strange and anti- 
strange quarks (m s = 0.45 GeV) are parameterized 
through a blast-wave with radial flow velocity v(r) = 
VQ-f/Ro. We compare the result for <fi mesons (m^ = 1.02 
GeV, = 0.05 GeV) calculated from resonance recom- 
bination of s and s, Eq. ([5]), with the direct blast- wave 
expression for the </> meson with identical flow field and 
temperature as for the quark input. Fig. [1] shows that 
both ways of computing the <j) spectra are in excellent 
agreement. We have reconfirmed that the results are in- 
sensitive to variations of the resonance width T when 
ensuring r < Q with Q = — (m s + rris) [l8j . 

To account for anisotropic flow fields which are suit- 
able for generating configurations reminiscent of hy- 
drodynamic simulations for fireballs in non-central (az- 
imuthally asymmetric) nuclear collisions, we adopt in 
this work the parameterization introduced by Retiere 
and Lisa (RL) [25j . It assumes longitudinal boost- 
invariance [26[ and parameterizes the transverse flow ra- 
pidity as a function of radius r and spatial azimuthal 
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where Tf is the freeze-out temperature at constant lon- 
gitudinal proper time 77, pi the chemical potential of 
particle i, gi the spin-isospin degeneracy factor, K\ a 
modified Bessel function and out- — pT/TfSiahp(r,(f> s ), 
(3t = mr /Tf cosh p(r,<fi s ). Equation ([9]) has been writ- 
ten in Boltzmann approximation which works well for 
sufficiently heavy particles; we will, however, use Bose 
distributions for pions. From Eq. Q we can calculate 
the elliptic flow of particle i as 



v l 2 ( PT ) 



Jo* # 



dNj 
P p T dp T d<j> p 



(10) 



FIG. 2: (Color online) Comparison of elliptic flow for <j>, D 
and J /if) mesons (top to bottom) computed from: (i) reso- 
nance recombination of strange quarks (symbols), and (ii) a 
direct application of the blast-wave expression at the meson 
level (lines). Both quark and meson distributions are based 
on the same anisotropic flow field constructed from the RL pa- 
rameterization with T c = 180 MeV, R x = 7.5 fm, R y = 8.5 fm, 
p = 0.66 and p 2 = 0.07. 



angle 4> s as 



p( r , <t>s) = r iPo + 92 cos(2</> b )] , 



(6) 



where po represents the average radial flow rapidity and 
P2 the anisotropy of transverse flow. The angle 4>b of the 
flow vector is generally tilted from the position angle <f) s 
(both angles are measured with respect to the reaction 
plane) which can be reflected by the relation 



tan (Ah = 



tan ( 



(7) 



This ansatz is motivated by the picture that the trans- 
verse boost is locally perpendicular to elliptical sub- 
shells 25]. The transverse flow rapidity, p(r,(f> s ), is as- 
sumed to increase linearly with the normalized "elliptic 
radius" , defined as 



r 2 sin 2 i 



s/Rl 



(8) 



We emphasize that this flow field incorporates space- 
momentum correlations that are considered a much more 
realistic representation of a collectively expanding source 
than the factorized ansatz of a "local" v 2 (independent of 
position) as often adopted in quark coalescence models 
calculations. 

Using the RL parameterization the differential momen- 
tum spectrum for a particle i directly emitted from the 
source takes the form 0, [2^, [27| (feed-down contributions 
are evaluated in App. [A} 



dN, 



p T dp T d4> P dy (27r) 3 

x j " rdr J d<t> s Kx(m T ,T,p T ) e * , (9) 



With the RL parameterization for an anisotropic flow 
field we are now in position to compare the elliptic flow 
generated for coalesced mesons in the RRM for locally 
equilibrated quark distributions, Eq. ([5J with the one 
directly obtained from applying the blast-wave expres- 
sion © at the meson level. We fix the parameters as 
quoted in the caption of Fig. [2] noting that at this point 
our aim is not to fit experimental data but to scruti- 
nize generic properties of the RRM. We emphasize again 
the nontrivial spatial dependence of the flow field figur- 
ing into the quark distributions in Eq. ([5]). To extend 
the scope of our comparison we consider, in addition 
to (j) mesons, also D {mo = 1.9 GeV, Yd = 0.1 GeV) 



and J/tp mesons (m 



3.1 GeV, T 



7/V 



0.1 GeV), as 



bound states of charm (m c = 1.5 GeV) and light quarks 
(m u ,d — 0.3, GeV) and of charm and anticharm quarks, 
respectively. Fig. [5] illustrates the comparison of the 
meson-w 2 (pT) for RRM and their direct blast- wave coun- 
terparts, showing again excellent agreement within nu- 
merical accuracy. It is quite remarkable that RRM gen- 
erates negative v-i at low px for the J/tp starting from 
strictly positive i>2(pt) for the underlying charm-quark 
distributions. A possibly negative i>2 is a well-known 
mass effect caused by the depletion of the heavy-particle 
yield at low momenta in the presence of sufficiently large 
collective flow Q- However, we are not aware of any 
quark coalescence model calculation that has reproduced 
this result. It demonstrates that the meson distributions 
resulting from RRM have reached local equilibrium and 
follow the collective flow of the source. We have checked 
that the negative v 2 is rather robust against variations of 
blast- wave parameters. By increasing po or p 2 the nega- 
tive value can be amplified. 

Combining the results from Figs. Q] and [U we conclude 
that hadron phase-space distributions obtained from 
resonance recombination precisely reflect the collective 
properties of a source in local equilibrium, explicitly doc- 
umenting the fact that the RRM achieves this through 
energy conservation and detailed balance. These insights 
also imply that the intermediate-pT regime, where v 2 sat- 
urates and shows an explicit dependence on quark num- 
ber, can not be in full equilibrium. 
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Group 


Tf 


Po 


P2 




Ry 


I 


180 


0.75 


0.072 


7.5 


8.9 


II 


110 


0.93 


0.055 


11.0 


12.4 



TABLE I: Parameter values for kinetic freeze-out tempera- 
ture, Tf in MeV, radial-flow rapidity po, rapidity-asymmetry, 
P2, and elliptic radii, J? x and R y in fm, for group-I and -II 
hadrons using the RL blast-wave expression for mid-central 
Au-Au collisions at RHIC. 



III. RECOMBINATION AND QUARK 
DISTRIBUTIONS AT HADRONIZATION 



Based on the verification that resonance recombina- 
tion maps the collective local-equilibrium properties of 
the quark source into hadron spectra, we can reverse the 
strategy and use this as a tool to extract quark distri- 
butions by a quantitative study of experimental data on 
hadrons which arise from quark coalescence. The gener- 
ality of this argument even allows to extend the analysis 
to baryons without having to explicitly calculate their 
spectra in RRM (which would be more involved than for 
mesons as it requires to recombine 3 quarks). As alluded 
to in the introduction, this procedure can only give access 
to quark distributions if the subsequent hadronic phase 
exerts a negligible distortion of the hadron spectra right 
after their formation around T c . We label such hadrons 
as belonging to "group I", and identify multi-strange 
hadrons (</>, 5, ft) as the prime candidates currently 
available from experiment. These hadrons have no well- 
established resonance excitations on the most abundant 
bulk particles (-7T, K and N) in the hadronic medium, 
while elastic i-channel exchange processes are suppressed 
by the OZI rule. One can therefore expect that multi- 
strange hadrons do not suffer significant rescattering in 
the hadronic phase, which is indeed supported by blast- 
wave fits to experimental data favoring kinetic decoupling 
not far from T c . The flow fields extracted for group-I 
hadrons therefore reflect the collective properties of the 
quark phase. 

We first determine the flow parameters in the RL 
parameterization by performing a blast-wave fit to the 
transverse- momentum spectra and elliptic flow of group-I 
hadrons to experimental data from STAR in mid-central 
(20-40%) Au-Au collisions. For definiteness (and easier 
comparison to the RRM calculations in Ref. [I(|), we 
fix the temperature to T c = 180 MeV; we do not utilize 
it as a parameter but obtain fits of satisfactory qual- 
ity with this value. The transverse-flow rapidity, po, is 
mostly governed by fits to the pr-spectra of <f>, 3 and f2 
hadrons, while the flow-asymmetry parameter, p2, and 
the fireball's spatial eccentricity, i? y /i? x , are driven by 
fitting the elliptic flow of S and <fi. The resulting RL 
blast- wave curves are shown with data in Fig. [3l and the 
central values of the parameters are collected in Tab. fl] 
As to be expected, the slopes of the spectra can be fit- 
ted well. The absolute normalizations are also repro- 
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FIG. 3: (Color online) Blast-wave fits to STAR data [lj, 
l2H - [3l| in semi-cental Au-Au( v / ijviv=200 GeV) collisions for 
transverse-momentum spectra (upper panel) and elliptic flow 
(middle panel) of 0, E + and (fl + fi)/2 using the RL ansatz 
with parameters specified in Tab. [I] The lower panel shows 
the resulting £l-to-(/> ratio on a linear scale, compared to STAR 
data [1<|. 



duced reasonably as well, which is confirmed by the Q- 
to-cj) ratio shown on a linear scale in the lower panel 
of Fig. [3J Degeneracy factors for spin and isospin are 
taken into account, but baryon chemical potential and 
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FIG. 4: (Color online) Upper panel: Equilibrated transverse 
momentum distributions of quarks on a hypersurface with 
temperature just above T c for a mid-central fireball at RHIC 
energies as predicted by the RRM model from group I data. 
Bottom panel: The same for the V2 of the equilibrated quarks. 



the strangeness suppression factor, 7 S , are neglected in 
these fits (they contribute at the 10% level which can 
easily be absorbed by fine-tuning the freeze-out temper- 
ature) . Note that the f2-to-0 ratio is quite different from 
linear for pr up to ^1.5 GeV, and that this ratio does 
not go to zero for px — ► as predicted by quark coales- 
cence models discussed in Ref . [H, H^] . This reiterates 
once more the importance of energy conservation and 
the thermal-equilibrium limit if on attempts to describe 
hadron spectra through quark coalescence in the low-px 
regime. 

We are now in a position to extract the quark distribu- 
tion functions at hadronization by assuming that: (i) the 
flow field for group-I hadrons is indeed the flow field for 
the hadronic phase in equilibrium just below T c (which 
is further supported by the fact that it is possible to use 
a freeze-out temperature of 180 MeV), and (ii) group-I 
hadron formation occurs via coalescence from a partonic 
phase (which is supported by the general CQNS charac- 
teristics of RHIC data). In Fig. 0] we not only include the 



extracted strange- and light-quark spectra and elliptic 
flow (which are directly involved in the RRM-based fits) 
but also those implied for the heavier charm and bottom 
quarks. Equilibration of the latter 2 quark flavors is cur- 
rently an open issue, but this may rather be a question of 
how far up in p t it applies. Qualitatively, the same limi- 
tation applies to light and strange quarks, which presum- 
ably enter the kinetic regime at pt — 1 GeV, where the 
data for the quark-number scaled v-i level off at a value 
of about 7-8%. Once available, D- and B-meson data 
will allow for similar estimates, and we provide here the 
underlying heavy-quark spectra as a reference for such 
an analysis. 

An extraction of strange- and light-quark spectra in a 
similar spirit, i.e., from data of what we refer to as group- 
I hadrons, has been attempted by in Ref. [32j based on 
a 1-dimensional instantaneous quark coalescence model. 
We reiterate here that such an extraction cannot be reli- 
ably applied at low px due to the simplifying assumptions 
inherent to such models, but it would certainly be inter- 
esting to make a comparison between both approaches. 



IV. SEQUENTIAL FREEZE-OUT AND BULK 
PARTICLES 

In URHICs the most abundant bulk particles, like pi- 
ons, kaons and protons (which we refer to as group-II par- 
ticles), undergo significant re-scattering in the hadronic 
phase. This is borne out of experimental blast-wave 
fits which require kinetic-freeze-out temperatures (e.g., 
T fo ~ 100 MeV for central Au-Au collisions at RHIC Q) 
which arc significantly lower than the chemical freeze- 
out temperatures (T c h — 160 — 180 MeV, extracted from 
observed ratios of different hadron species). Note that 
the difference in the kinetic and chemical freeze-out tem- 
peratures translates into a factor 5-10 different energy 
densities [Hj]. Theoretically, the sequential freeze-out 
systematics can be readily understood from large reso- 
nant cross sections operative for group-II particles in the 
hadronic phase, e.g., 7T7t — > p — > tttt, ttK — > K* — > ttK 
or ttN — > A — > irN. The values of these elastic cross sec- 
tions reach up to 100-200 mb and are thus 1-2 orders of 
magnitude larger than chemistry-changing inelastic ones 
(e.g., 7T7r — )• KK). Further direct evidence for an ex- 
tended duration of hadronic reinteractions follows from 
the well-established low-mass dilepton enhancement in 
URHICs whose magnitude and spectral shape requires 
several generations of p mesons to form with a largely 
broadened spectral shape due to a strong coupling to the 
medium [34j . 

In order to include group-II particles in an analysis of 
the thermal regime of anisotropic flow fields there is thus 
no choice but to introduce an additional source parame- 
terization at a much lower temperature than for group-I 
particles. While this generally detaches group-II parti- 
cles from coalescence mechanisms at hadronization, we 
remark that their production through resonance recom- 
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FIG. 5: (Color online) Blast- wave fits using the RL fireball pa- 
rameterization to proton and kaon pr-spectra (upper panel) 
and elliptic flow (lower panel) in mid-central Au-Au collisions 
at RHIC. The flow- field source parameters are listed in Tab.|l] 
Experimental data are taken from Refs. [3. l35l. l36l] . 



bination in equilibrium at T c by no means contradicts 
but rather facilitates local equilibrium until a much lower 
bulk freeze-out. Group-II hadrons at low px are there- 
fore not suitable to infer direct statements about quark 
recombination. Our aim rather is to test whether the 
minimal assumption of a sequential freeze-out is suffi- 
cient to produce the general scaling properties observed 
in the data, and, if so, what the key features are and 
how sensitive this picture is to variations of the freeze- 
out configuration. We note in passing that the situation 
is presumably quite different at intermediate pr where 
hadronic re-scattering is insufficient to maintain equili- 
bration and one is rather dealing with the (transition to) 
a kinetic regime. 

Toward this end we again employ the RL blast-wave 
ansatz to parameterize an anisotropic flow field which 
results in a good fit to the px spectra and elliptic flows 
of group-II particles. As in the group-I case, we decide 
to fix the freeze-out temperature at a value compatible 
with earlier blast- wave extractions, at Tf Q = IIOMeV. 
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FIG. 6: (Color online) Parameter contours for p2 and i? y /i? x 
enclosing the regions of consistency with RL blast-wave fits 
to pt spectra and V2 in the thermal regime for both freeze- 
out groups. The temperature, T, and average radial-flow ra- 
pidity, po, have been fixed at (T, po) = (180 MeV, 0.75) and 
(110 MeV, 0.93) for group-I and -II, respectively. 

Then we determine po, p2 and the eccentricity of the 
fireball from fits to the measured spectra and V2 of pro- 
tons and kaons for the same centrality selection of Au- 
Au collisions. The results are displayed in Fig. [5j cor- 
responding to the parameter values listed in Tab. |U At 
kinetic freeze-out, effective chemical potentials for stable 
particles (pion, kaon, nucleon) are required to corr ectly 
predict the absolute normalization of their yields [33| . 
Neglecting chemical potentials therefore implies that we 
can only fit the shape of the spectra which is sufficient 
for our purposes here. In Boltzmann approximation, the 
chemical potential appears as a constant overall fugac- 
ity factor which does not affect the V2 either. The pions 
are treated with Bose statistics which entails an appre- 
ciable (moderate) enhancement of the px spectra (^2) at 
Pt ;$ 0.5 GeV. At the same time, feed-down contribu- 
tions are also significant in this regime. A more detailed 
discussion of the latter is given in App. [5] 

It is instructive to elaborate on the qualitative influ- 
ence of the fit parameters on the observables. As is well- 
known from more simple blast- wave fits, the interplay of 
radial flow (po) and temperature (T) can be disentan- 
gled by fits to hadrons with different masses. With these 
2 parameters fixed, an increase in P2 leads to a rapid in- 
crease in V2(pt) at relatively large pr but to rather small 
changes at low px- On the contrary, increasing R y /R x 
results in a more uniform increase of V2(pt) over a wide 
range of px- As a result, a relatively large P2, combined 
with a relatively small eccentricity, i.e., R y /R x close to 
unity, are favored in fitting the data. The rather tight 
correlation between those two parameters can be exhib- 
ited in pertinent contour plots indicating the "allowed re- 
gions" for obtaining an acceptable fit to data given their 
error bars, see Fig. [5] The central values of the contours 
for the two particle groups correspond to the values listed 
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in Tab. HI The contour for group-II hadrons is tighter due 
to the larger number of data points for bulk particles. We 
first note that the contours clearly exclude each other 
corroborating the significant evolution that the collec- 
tive flow field is undergoing in the hadronic phase. More 
specifically, the anisotropy parameter R y /R x for group-II 
particles is smaller than for group-I, which is in line with 
the general expectation that the fireball becomes more 
spherical in the expansion. However, the flow anisotropy 
P2 extracted from the RL blast-wave fits favors a decrease 
from hadronization (group-I) to kinetic decoupling of the 
bulk (group-II) . This is contrary to the expectation from 
hydrodynamic expansion that pi — (p£ — Py)/2 should 
increase if the in-plane acceleration is always larger than 
the out-of-plane acceleration (p£ and p y are the surface- 
flow rapidities along the x- and y-direction, respectively; 
recall Eq. ([6])). We cannot resolve this question within 
the present blast-wave model. A full hydrodynamic sim- 
ulation could give more insight, possibly requiring sub- 
stantial viscosity effects in the later stages of the hadronic 
evolution. 

To summarize this and the previous section, we 
have produced a space-dependent parameterization of 
anisotropic flow fields which, when implemented into a 
rather simple thermal blast-wave description of locally 
equilibrated matter, provides a good overall fit to a large 
body of hadron data in semi-central Au-Au data in the 
low-pT regime. As a key ingredient we have introduced 
the well established concept of sequential freeze-out in- 
cluding just two groups but with clearly distinct flow 
fields. 



V. KE T SCALING OF v 2 

A signature result of ideal hydrodynamical simulations 
of URHICs is the mass splitting of the hadron V2 when 
plotted as a function of hadron pt, with a larger V2 for 
the lighter particles at a given px- On the other hand, 
if the hydrodynamic results for the hadron vi are plot- 
ted versus KEt (at a given freeze-out temperature), the 
splitting nearly vanishes, but not entirely. In fact, a small 
reverse mass ordering of V2 is found instead, as pointed 
out in Refs. [H,!!?]]. In other words, the hydro-induced 
mass splitting of «2 (pt) between light and heavy particles 
is not large enough to render their V2{KEt) curves degen- 
erate. The sequential freeze-out picture advocated above 
may hold an explanation to this apparent discrepancy 
between hydrodynamic freeze-out (at fixed temperature) 
and experimental data (which exhibit "perfect" KEt 
scaling): Random thermal motion tends to isotropize pt 
distributions and thus counteracts ( "dilutes" ) the collec- 
tive motion in producing nonzero t>a. Thus one should 
expect that a gap in freeze-out temperatures between, 
say, (light) pions and (heavy) Q baryons increases the 
mass splitting in a V2(pt) plot (the thermal motion of 
pions is suppressed at lower temperatures), just enough 
to produce degeneracy in a V2(KEt) representation. This 
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FIG. 7: (Color online) KEt- and n ? -scaling of V2 for differ- 
ent hadrons obtained within the sequential freeze-out scenario 
with blast-wave parameters given in Table [I] The result of an 
slightly later freeze-out for pions at temperature 100 MeV is 
also shown. 

appears to be a conspiracy of circumstances that leads to 
the empirically observed i£&r-scaling. However, the un- 
derlying mechanism (largely different hadronic rescatter- 
ing cross sections) is general in the sense that it applies to 
different centralitics, system sizes and even collision en- 
ergies, presumably coupled with a rather general (local) 
freeze-out criterion. It therefore could be a candidate for 
the origin of the widely observed KEt scaling. In the 
following, we quantitatively test this idea by analyzing 
the scaling properties within our sequential freezeout 2- 
source ansatz. 

With the RL flow-field parameterizations and freeze- 
out temperatures for the two sources of group-I and 
group-II particles summarized in Tab. U we proceed to 
display the rt g -scaled hadron-t^ as a function of the 
n g -scaled hadron- KEt in Fig. [7] We find that all the 
hadrons-«2 curves follow a generic scaling very well (with 
slights deviations for the pions, on which we comment be- 
low), very similar to experiment. This agreement is, of 
course, not unexpected given that our results follow from 
fits to data, together with the notion that the data satisfy 
the scaling in the first place (or even led to its discovery) . 
The nontrivial insight here is that our "minimalistic" ap- 
proach of a sequential freeze-out scenario coupled with a 
rather simple parameterization of an anisotropic collec- 
tive flow field, is sufficient to establish the scaling. We 
also recall that n g -scaling in the equilibrated low-py re- 
gion is not a direct consequence of quark recombination, 
at least not for group-II particles. To better exhibit the 
significance of the sequential freeze-out we compare in 
Fig.[5]the generic scaling curve extracted from Fig.[7]with 
V2 results where the (rather heavy) f2 baryon (group-I) 
is frozen out from the source at Tf = 110 MeV (as some- 
times done in hydrodynamic calculations with a single 
freeze-out), and the (rather light) pions and kaons are 
frozen out from the source at T c = 180 MeV (as advo- 
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FIG. 8: (Color online) Comparison of the generic KEt- 
scaling behavior of hadron-i>2 extracted from Fig. [7] (thick 
solid line) to results for interchanging the freeze-out configu- 
rations for group-I and group-II particles. Specifically, the Q 
baryon (dash-dotted line) is frozen out from the late source at 
Tf = 110 MeV, while pions (dotted line) and kaons (dashed 
line) are frozen out from the early source at T c = 180 MeV 
(source parameters as quoted in Tab. [I]). 



cated in recombination models that neglect the hadronic 
phase) . One finds a considerable breaking of both KEt- 
and n g -scaling. For the heavy particles, the increase in 
collective flow at the low freeze-out temperature results 
in an appreciable depletion of their at low KEt (mass 
ordering effect), thus spoiling their compatibility with 
scaling. On the other hand, for pion and kaon freeze- 
out at higher temperature the increased thermal motion 
in the local rest frame dilutes the collective elliptic flow 
considerably and suppresses their vi below the scaling 
curve especially toward higher KEt- 

Let us return to the deviations of the pion V2 from the 
scaling curve which become noticeable for KEr/n q > 
0.3 GeV in Fig. [7] One may wonder whether feed-down 
corrections from heavier particles (p — > 7T7t, uj — > Sir, 
A — > Nir, etc.) could improve the situation. This ques- 
tion is addressed in detail in Appendix[Al with the finding 
that the combined effect of the decay pions increases the 
pion v-i at pr > 0.5 GeV only a little (by a few percent). 
It turns out, however, that a slight decrease of the pion 
freeze-out temperature from 110 MeV to, e.g., 100 MeV 
increases the pion V2 preferentially at higher KEt, re- 
sulting in good agreement with the general scaling curve. 
In doing so we have not performed a possibly small re- 
adjustment of the flow parameters to attempt a more 
quantitative conclusion. The key mechanism is again a 
further suppression of the random thermal motion of the 
pions, thus augmenting the relative importance of the 
collective motion. A 5-10 MeV smaller freeze-out tem- 
perature is certainly compatible with experimental data, 
and also not unreasonable from the theoretical point of 
view. Pions have large resonant cross sections on all ma- 



jor particle types of the bulk (pion, kaon and nucleon), 
and their mass being the smallest makes it plausible that 
their momentum spectra are frozen the latest (mass or- 
dering of thermal relaxation times) . 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have analyzed the scaling properties 
of hadron elliptic flow in ultrarelativistic heavy-ion col- 
lisions, which have received considerable attention since 
the experimental data suggest the existence of a univer- 
sal behavior possibly related to a collectively expanding 
partonic source. We have focused our investigations on 
the low-pr regime (pr < 2 GeV) where RHIC data (and 
hydrodynamic analysis thereof) support the assumption 
of a locally equilibrated medium. 

In the first part of this work we have pursued the ear- 
lier idea that hadronization of the putative partonic sys- 
tem at RHIC proceeds through coalescence of quarks and 
antiquarks. Its implementation in the thermal regime re- 
quires a recombination model which obeys 4-momentum 
conservation and detailed balance to encode the correct 
equilibrium limit. The Resonance Recombination Model 
(RRM) satisfies these requirements. The flow field has 
been modeled by a blast-wave ansatz with azimuthal 
asymmetries which give rise to space-momentum corre- 
lations characteristic for a collectively expanding source. 
Using RRM we have shown that an expanding parton 
source converts into hadron spectra with the same collec- 
tive source properties, which, to our knowledge, has not 
been achieved in existing models of quark coalescence. 
We have utilized this connection to extract quark dis- 
tributions just above the critical temperature for Au-Au 
collisions at RHIC from measured pr spectra and elliptic 
flow of multi-strange hadrons which are believed to un- 
dergo negligible final-state interactions in the hadronic 
phase (group-I particles). 

Bulk particles with large hadronic cross sections (e.g., 
pion, kaons and nucleons) cannot be approximated by 
decoupling at the hadronization transition. We have in- 
cluded these particles (group-II) in our analysis by con- 
structing another flow field at the typical kinetic freeze- 
out temperature of 110 MeV, and fitting the blast-wave 
parameters to their pt spectra and elliptic flow. Group-II 
hadrons are therefore not really sensitive to quark recom- 
bination. However, we find that the sequential freeze- 
out picture with just two sources is well compatible with 
the empirically observed kinetic-energy and constituent- 
quark number scaling (where the latter is not directly 
linked to recombination) . The sequential freeze-out natu- 
rally corrects for a slight breaking of KEt scaling present 
in a hydrodynamic flow description for a uniform freeze- 
out. In particular, a smaller temperature is instrumental 
in suppressing the random thermal motion of the light 
pions and in this way bring up their V2 toward the uni- 
versal curve. In the thermal regime, CQNS of elliptic 
flow simply follows from the universal linear curve re- 
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suiting from KEt scaling (i.e., it is not directly linked to 
recombination) . 

An important task of future work will be an exten- 
sion of our studies to the intermediate-pT region. Here, 
the onset of the kinetic regime occurs and thus off- 
equilibrium effects become important. While this in- 
creases the sensitivity to quark degrees of freedom, it also 
requires a more detailed knowledge of the reinteractions 
of particles, both in the QGP and hadronic phase. Initial 
studies using Langevin simulations for strange and heavy 
quarks have been reported [20l| . but the scaling proper- 
ties found there clearly require a deeper understanding, 
including the effects of baryon formation and a more mi- 
croscopic description of the elastic (as well as radiative) 
interactions in the QGP and hadronic matter. 
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Appendix A: Feed-Down to Pions From Resonance 
Decays 

Hadron-resonance decays are known to increase pion 
yields significantly for transverse momenta up to pr — 
0.5-1 GeV [38]. Therefore, it is mandatory to evaluate 
the effect of the corresponding feed-down not only on the 
Pt spectra but also on the elliptic flow of the pions [H, 
|40| . In particular, we want to understand how feed-down 
affects the pion v-i in the context of the iffir-scaKng for 
pion freeze-out at a temperature of Tf a = 110 MeV where 
some deviation from protons and kaons in group I has 
been identified (recall Fig. [7]). The total differential pion 
spectrum can be decomposed as 

dK ot = dN^ | rfAW 
p T dp T d(j) p dy p T dp T d4> p dy ^ p T dp T d<f> p dy ' 

where the summation is over all resonances R with sig- 
nificant branching fraction into final states containing pi- 
ons. In our decay simulation, we include p, oj, K*, A, r\ 
and Kq resonances, whose anisotropic momentum spec- 
tra, together with that of direct pions, are obtained from 
the RL blast-wave calculation with the same parameters 
as those for group II in Tab. [I] The resonance decays 
which contribute to the final spectra occur at or after 
freeze-out. However, the absolute number of the feed- 
down contribution, which is important to determine the 
relative contribution to the total pion-U2, depends on the 



chemical potentials of the resonances at kinetic freeze- 
out. In the interacting hadronic phase effective chemical 
potential need to be introduced for hadrons stable under 
strong interactions (pion, kaon, t], nucleon, etc.), as to 
conserve their total numbers which are frozen at chemi- 
cal freeze-out (4lTl43l |. Part of this procedure is to keep 
the strong resonances in relative chemical equilibrium, 
e.g., \ip = 2/ijr or = + Mtt- At RHIC, the con- 
servation of antibaryon number is of particular impor- 
tance in the overall chemistry of the hadronic phase, and 
thus for the actual values of the effective chemical po- 
tentials [§1, |H, |44T - |46T | . The values of the chemical poten- 
tials of the resonances at Tf Q are quoted in the caption of 
Fig. El they are taken from Ref. |33[ where they are com- 
puted from an iscntropic evolution starting at chemical 
freeze-out at T — 180 MeV (consistent with our param- 
eters for group- 1 freeze-out). Note that the Kq decays 
through weak interactions but carries its own chemical 
potential (as does the rj). 

The p^-spectra of pions from various resonance decays 
are compiled in Fig. [S] The decay pions generally have 
steeper p^-spectra than directly produced pions, except 
for pions from p decays, similar to the observations made 
in Refs. [39|, H3|. This is essentially a phase-space effect: 
pions from resonances with a small Q-value (the differ- 
ence of the resonance mass and the sum over masses of its 
constituent quarks) tend to have a steeper px spectrum. 
The phase-space also affects the v-i of decay pions. The 
momenta of pions from resonances with large Q values 
(like the p) tend to dilute the V2 more due to the ran- 
dom orientation of the decay momentum, most notably 
at low momenta (where the decay momentum in the par- 
ent's rest frame is comparable to the parent momentum). 
On the other hand, pions from resonances with small Q 
values (like the A) are essentially aligned with the par- 
ent and thus preserve the elliptic flow of the parent. For 
those pions, the V2 could be larger than that of the parent 
at the same pt, since the momentum of the resonance is 
shared by several daughters. However, since the dom- 
inant source of secondary pions are p decays, the total 
pion-^2 might not be much affected by resonance decays. 

The comparison of the total pion-i>2 to the direct-pion- 
V2 confirms this expectation, see Fig.JTU] Therefore, feed- 
down contributions are not sufficient to resolve the dis- 
crepancy between the direct-pion v-i and the generic scal- 
ing curve in the ifE-r-dependence of V2- In comparison 
to previous estimates we note that in Ref. [39| a larger 
enhancement of the total pion-i>2 over the direct-pion i>2 
was reported. In that work, the resonance decays were 
evaluated at T c ; subsequent hadronic rescattering was ne- 
glected and the direct pions were computed from a quark 
coalescence model which does not conserve energy. Both 
of these approximations may be unreliable for pions at 
small pt, due to their large interaction cross section and 
binding energy, respectively. The influence of resonance 
decays on pion elliptic flow was also studied in Ref. [3] , 
where the authors assumed that direct pions follow the 
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FIG. 9: (Color online) Transverse momentum spectra of pions 
from various resonance decays compared to direct pions (solid 
line) at Tf = 110 MeV. The chemical potentials for the par- 
ent particle are resonances are /i p = 180 MeV, fiA = 460 MeV, 
fMK* = 250 MeV, n u = 270 MeV, n„ = 170 MeV, result- 
ing from a partial chemical equilibrium scenario in which 
chemical potentials for stable particles are fj,^ = 90 MeV, 
UK = 160 MeV and fx N = 370 MeV. 



FIG. 10: (Color online) Total pion-«2 (dashed line) compared 
to direct pion V2 (solid line) at T{ = 0.11 GeV. 



KEx-scalmg and then found the scaling still holds at the 
10% level when resonances are included. 
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